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Fiiud  Technical  Report  on 
Numerical  Simulation  of  the  Collapse 
of  an  Underwater  Explosion  Bubble 

OMR  Contract  N00014-87-K-0428 

> 

Puroose:  The  goal  of  ^s  project  waa^demonstrate^the  usefulness  of  front  tracking.  (J. 
GUmm.  eL  al.,  Courant  Institute,  NYU]hts  a  tool  for  studying  the  fully  compressible,  axially  sym¬ 
metric  expansion  and  collapse  of  an  underwater  explosion  bubble  in  two  spatial  dimensions.  In 
the  pilot  project  supported  by  this  contract,  a  one  dimensitmal  version  of  the  ftont  tracking  code 
was  implemented  to  compute  statically  symmetric  oscillations^ 

Final  status  and  conclusions:  SA  one  dimensioiutl,  tracked  random  choice  code  [1]  was 
adapted  for  this  pilot  project,  and  validated  against  published  computational  and  experimental 
results.  Specifically,  comparisons  were  made  with  untracked,  random  choice  calculations,  by 
Saito  and  Glass  [2]  and  Botes  and  Holt  [3]  and  by  calculations  of  Mader  [4]  using  the  code  SIN. 
These  validation  studies  were  ptesetued  in  the  Appendix  to  the  previously  submitted  Annual 
Report  for  1987  for  this  contract  A  copy  of  this  annual  report  is  enclosed. 

These  validation  studies  identified  several  improvements  requited  by  the  tracked  algorithm. 
These  included  1)  the  necessity  for  realistic  equations  of  state  for  the  explosion  products  (stif¬ 
fened  gamma  law  gas  equations  of  state  where  used  in  the  study),  2)  the  implementation  of  effec¬ 
tive  boundary  condititms  that  would  take  into  account  the  progress  of  the  primary  shock  wave 
into  the  water  without  the  necesrity  of  keeping  it  within  the  computational  domain,  3)  improved 
treatment  of  the  divergence  at  the  bubble  center,  and  4)  correa  initiation  (i.e.  velocity  profile)  of 
the  explosion  remairts.  — 
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Plans  were  fonnulated  for  overcoming  each  of  these  problem  areas,  and  were  discussed  in 
the  aforementioned  annual  report  Briefly  they  were:  1)  implementation  of  the  HOM  equations 
of  state,  which  were  to  be  input  to  a  Reimann  problem  solver  developed  for  general  tabulated 
equation  of  state  data  by  J.  Scheuerman  at  tlw  COurant  Institute,  2)  using  the  ideas  of  T.  Hagstrom 
[S]  and  others  to  achieve  effective  far  field  boundary  coiditions  once  the  primary  shock  leaves 
the  domain  of  computation.  3)  investigating  the  analytical  treatment  of  Noh  as  a  possible  method 
for  handling  the  divergences  at  the  bubble  origia  4)  implementation  of  a  model  for  the  reactive 
phase  of  the  explosion  in  one  of  the  marmers  discussed  by  Mader  [4], 

This  project  was  terminated  at  the  end  of  the  pilot  projea,  when  further  research  funding 
was  discontinued.  Other  than  the  Appendix  in  the  previously  submitted  annual  report,  no  techni¬ 
cal  report  or  publication  resulted  from  this  project. 
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